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ABSTRACT 

Recent work suggests that Type la supernovae (SNe) are composed of two distinct populations: prompt and 
delayed. By explicitly incorporating properties of host galaxies, it may be possible to target and eliminate 
systematic differences between these two putative populations. However, any resulting /^oif-calibration shift 
in luminosity between the components will cause a redshift-dependent systematic shift in the Hubble diagram. 
Utilizing an existing sample of 192 SNe la, we find that the average luminosity difference between prompt 
and delayed SNe is constrained to be (4.5 ± 8.9)%. If the absolute difference between the two populations 
is 0.025 mag, and this is ignored when fitting for cosmological parameters, then the dark energy equation of 
state (EOS) determined from a sample of 2300 SNe la is biased at ~ la. By incorporating the possibility of 
a two-population systematic, this bias can be eliminated. However, assuming no prior on the strength of the 
two-population effect, the uncertainty in the best-fit EOS is increased by a factor of 2.5, when compared to 
the equivalent sample with no underlying two-population systematic. To avoid introducing a bias in the EOS 
parameters, or significantly degrading the measurement accuracy, it is necessary to control the post-calibration 
luminosity difference between prompt and delayed SN populations to better than 0.025 mag. 
Subject headings: cosmology: observations — cosmology: cosmological parameters — supernovae: general 
— surveys 



1. INTRODUCTION 

The discovery of the accelerating expansio n of the uni- 
verse dRiess et al]|1998t iPerlmutter et al.ll 19991) has led to an 
explosion of interest in the underlying physics responsible for 
this acceleration. A favored model characterizes the accel- 
eration by an unknown energy density component, dubbed 
the dark energy. While there exist a variety of probes to ex- 
plore the nature of this dark energy, one of the most com- 
pelling entails the use of type la supernovae (SNe hence- 
forth) to map the expansion history of the universe. Sev- 
eral present and future SN surveys are aimed at constrain- 
ing the dark energy equation-of-state (EOS) to better than 
10%. With increasing sample sizes, SN distances can poten- 
tially provide multiple independent estimates of the EOS of 
dark energy when binned in reds hift ( Huterer & Coora vl2005l : 
ISullivan. Cooray. & Hold 120071: ISarkar et alJl2008al) . Given 
the importance of dark energy measurements, it is then use- 



ful to quantify various systematics that impact SN cosmology 
dHui & Greendl200aiCooray & Caldwelll u2006HCc~ 
I20C 



.oorav 



et alJ 



2006t ISarkar et al.ll2008bl) . 

Recently, suggestions have been made that the SN popu- 
lation consists of two components, with a "prompt" compo- 
nent proportional to the instantaneous host galaxy star for- 
mation rate, and a "delayed" (o r "extended") component 
that i s delayed by several Gyrs dHamuv et all [l995l iLiviol 



2000t IScannapieco and Bilds enll2005t iMannucci et al 
Sullivan et al. H2006tlStrovinl 2007). The former is expected 



to be more luminous, and thus, prompt SN lightcurves are 
broader than those of the delayed popul a tion. B y classifying 
SNe by host galaxy type, iHowell et all (120071) found a pre- 
calibration intrinsic luminosity difference of ~ (12 ±4)% be- 
tween the two components, based on a difference of (8.1 ± 
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2.7)% in the width of lightcurves. Since the SN lightcurves 



are used to calibrate the intrinsic lumino s ity JPhi llips 



Riess et all 



Prieto et all 



1996; Perlmutteret a 



11997b iTonrvetalJ 



1993 



2003 



2006; iGuv et al.ll2007l;ljha et alj|2007l) . a system 



atic difference in intrinsic luminosity could conceivably be 
calibrated out, if the SN lightcurve calibration relation is the 
same for both populations. 

However, it is unclear whether the full intrinsic difference 
in luminosity between the two populations is captured by a 
calibratable difference in the lightcurves. A residual in the 
calibrated luminosity could potentially remain, leading to a 
redshift-dependent shift in the Hubble diagram, and system- 
atic errors in the best-fit cosmological parameters. For exam- 
ple, it is likely that the intrinsic colors of Typ e la SNe are 
not u niquely determined by light-curve shape dConlev et all 
120081) . Even if this is not the case, differences in intrin- 
sic color between the populations might introduce system- 
atic differences in post-calibration luminosity (e.g., through 
differences in the extinction corrections). We model the 
two-population systematic, constraining the magnitude of the 
effect with current data. With large SN samples it may 
be possible to estimate the magnitude of the systematic di- 
rectly from the data (e.g., by correlating observ ed SN bright- 
nesses with properties of the host galaxies (lHamuy et alJ 



nesses with prope rties or the host g alaxie s (lHamuy et alJ 
19965 lRiessetal.ll 19991: ISullivan et al.ll2006llJha et al.ll2007h 
Gallagher et alJ 120081) ). We quantify the level of calibration 



required to avoid significantly degrading the determination of 
the dark energy EOS. 

The paper is organized as follows: in §2 we discuss a model 
for incorporating a two-population systematic residual lumi- 
nosity into the Hubble diagram. In §3 we investigate the possi- 
bility of detecting this systematic from both current and future 
SN data, and the impact on dark energy parameter estimation. 

2. TWO SN POPULATIONS AND THE HUBBLE DIAGRAM 

The use of SNe as standardizable candles to constrain the 
dark energy EOS is based on the fundamental assumption that 
the lightcurves of all individual SNe can be calibrated. The 
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FIG. 1 . — Left panel: Correlation between So and w, corresponding to a wCDM model, for the 192 SN dataset discussed in the text. There is a strong degeneracy 
between the two-population systematic, So, and the dark energy EOS, w. This leads to a potential bias in the measurement of the EOS, and increased errors when 
the fits are marginalized over So- Right panel: Histograms showing the distributions of the best-fit w, from 200 mock data sets with an inherent So = 0.025, 
after marginalizing over the WMAP wCDM 5-Year priors. The shaded histogram represents the case where So is allowed to vary freely, while for the hatched 
histogram the two-population systematic is ignored (So = 0). The dot-dashed line corresponds to the case where the data is fit assuming a So = 0.025 ± 0.025 
prior. For clarity we omit the underlying histogram. 



lightcurve shape-intrinsic luminosity relation is derived us- 
ing low-z samples of SNe, and it is assumed that this relation 
is applicable to the higher- j population. However, if the SN 
population is non-uniform wit h redshift, the potential for sys- 
tematics must be considered dHowellll200lHMannucci et alj 
l2006tlScannapieco and Bildstenll2005HSullivan et al.ll2006l) . 

According to the two-p opulation model of 
IScannapieco and Bildstenl (120051) (henceforth SB), the 
prompt SNe track the instantaneous star formation rate, 
M*(f), and dominate the total SN rate at early times (high 
redshifts) when star formation is more active. The rate of 
the delayed component scales proportionally to the total 
stellar mass at a given instant, M*(f), and thus delayed SNe 
dominate at late times (low redshifts). The total SN rate can 
be written as 



(lOOyr)" 1 



M»(0 



10 10 M, 



(3 



+ B 
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lO ll) M Gyr" 
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where A and B are dimensionless constant s. We take the 
star f ormation rate proportional to e"'A 2Gv '') dMannucci et al.l 
120061) . This model is likely to be a simple approxima- 
tion t o a true, smooth underlying distribu t ion of delay 
times dTotani et al.ll2008t iGreggio et al.1 120081: iPritchet et al.l 
120081) . For simplicity we restrict ourselves to the two- 
population model, although in j j3.2| we discuss sensitivity to 
this assumption. 

We allow for a redshift-independent residual difference 
(AL) in the calibrated absolute luminosities of the prompt 
(L P ) and delayed (L D ) SNe: L P = L D + AL. We take the de- 
layed and prompt fractions of the total SN population to be 
f D (z) and f P (z), respectively (f D (z)+fp(z) = 1). 

The average of the distance moduli of all the SNe in a given 
redshift bin can be written as 



((m-M)} =51o gl 



(—)■ 
VMpc / 



+25-(2.51og 10 



-vef 



(2) 



where L re f is the reference luminosity corresponding to the 
reference absolute magnitude M re r given by M re r = /dMd + 
fpMp, where f t is the redshift-average of fi(z) over the low 
redshift range where the calibration is done. The residual 
systematic correction to ((m-M)} can be expressed as ((am — 



M)) res = S D f D (z)-2.5log m (Lp/Lref), where S D = 1.086 ln(l + 
AL/Lo) is the two-population bias. Note that one can also 
express ((m-M)) les relative to the prompt component, using 
Sp instead of 5d, and fp instead of fo- This involves an overall 
sign change, but the final results are unaltered. 

Using the x 2 -statistic, we fit the Hubble diagram with a 
modified form of the distance modulus which includes a pos- 
sible residual: 



(m-M) fid (z) = 51og 



10 



ddz) 
Mpc 



+ 25 + M + S D f D (z). (3) 



We have absorbed the redshift-independent term in ((m — 
M))ms into the "nuisance parameter", M.. Note that Hq is 
incorporated into c/ £ instead of M.. The redshift-dependent 
factor foiz) can be determined based on our knowledge of the 
star-formation history, and So must be estimated directly from 
the SN data. This residual systematic must now be marginal- 
ized over, and will affect cosmological parameter estimates. 

3. RESIDUAL SYSTEMATIC AND PARAMETER ESTIMATION 

3.1. Existing data 

To study the extent to which existing data may be affected 
by a residual systematic in the luminosity difference between 



set of 192 SNe ( 


Davis et al.ll2007t IWood-Vasev et all 12007: 


Riess et al. 


|2007 


; Astier et al.1 120061 lHamuv et al.1 1996a; 


Riess et al. 


19991 


Jha et alj2006l). We also include two barvon 



acoustic o scillation (BAO) dis tance estimates at z = 0.2 and 
z = 0.35 dPercival et al.l 120071) . and the dimensionless dis- 
tance to the surface of last scattering R = 1.710 ±0.019 
dKomatsu et al.l 120081) . We take a flat ACDM cosmologi- 
cal model and marginalize over M., with W MAP priors of 
H = 71.9 ±2.6 and n„,h 2 = 0.1326 ± 0.0063 dKomatsu et alj 
12008b . These priors are independent of the SN data. For sim- 
plicity we do not incorporate the correlation between £l m h 2 
and Ho. Such a correlation will not qualitatively alter our re- 
sults, but will need to be taken into account once precision 
data becomes available. 

The distribution function for So from existing data has a 
positive mean value (So = 0.049), implying that the post- 
calibration luminosity of the delayed population is dimmer 
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FIG. 2 . — Left panel: The fraction of the delayed component, /b, as a function of redshift. The solid line corresponds to the fiducial A and B values 
from SB I Scannapieco and Bildsten 2005) that are used to generate the mock catalogs. The dashed and dot-dashed lines correspond to {A — AA, B + AB} and 
\A+ AA,B— AB}, re spectively, where AA and AB are the lcr uncertainties on A and B from SB. The triple-dot-dashed line depicts the recent estimate of fo from 
lAubourg et al. 1 2007J). Right panel: Histograms showing the best-fit values for the dark energy EOS parameter, w, for 200 mock SN datasets with So = 0.025. 
The solid curve shows the Gaussian fit to P(w) for the case where we use the same fb(z) to generate the mock catalogs and fit the data. The other Gaussian 
curves correspond to the best-fit values of w when we fit the mock data with the wrong functional form for /d(z). The mocks are generated with the fiducial 
model, and the line-types represent the model used in the fit, in accordance with the left panel. For clarity we only show two underlying histograms, for the cases 
{A — AA ,B + AB} (shaded) and {A + AA,B— AB} (hatched). A ~ 30% shift in fo (left panel) corresponds to a percent-level bias in w (right panel). The extreme 
model of Aubourg et al. 12007) corresponds to an almost redshift-independent fo (see left panel), and thus the bias is similar to that in the So = case. 



than the prompt population by ^5%. The standard deviation 
of this 5d is a = 0.097, which suggests that current SN data 
is consistent with th e absence of a two-population bias. The 
iHowell et all d2007l) result of - (12 ± 4)% difference in the 
intrinsic luminosity is a pre-calibration difference. We find a 
~ (5 ± 9)% difference in the posf-calibration luminosity be- 
tween the two types of SNe. While this uncertainty is larger 
than the pre-calibration value, it is estimated directly from the 
Hubble diagram, and is independent of the empirical stretch- 
luminosity relation. This important consistency check will 
improve as the SN sample sizes increase. 

To study the impact of a potential luminosity differ- 
ence on the measurement of the dark energy EOS, we 
also fit the same data (SN+BAO+CMB) assuming a wCDM 
model, and take the corresponding WM AP+HST priors Hp = 
72.1 ±7.5 and VL m h 2 = 0. 1329 ± 0.0066 dKomatsu et al.ll2008l: 
iFreedman et al]l2~001l) . We first consider two extreme cases: 
ignoring the two-population systematic (6d = 0), and allow- 
ing for a completely unconstrained systematic (So with no 
priors). With So free, our analysis yields a best-fit time- 
independent EOS parameter w = -0.986 ±0.1 80, with a best- 
fit 6 D = 0.040 ±0.279. If we set S D = 0, the same data pro- 
vides a constraint of w = -0.960 ± 0.066. By incorporating a 
two-population effect in the fit, the errors on the best-fit EOS 
degrade by a factor of ~ 3. In addition, we note a shift in the 
best-fit value of w between the two cases. The evolution in 
the ratio between prompt and delayed SNe, as a function of 
redshift, can mimic dark energy. If there is a residual differ- 
ence in the calibrated luminosity, ignoring 5p might bias the 
dark energy estimate. As we show in the left panel of Fig- 
ure 1, the increase in the error of the best-fit EOS is due to the 
degeneracy between dark energy (w) and the two-population 
systematic (8p). 

Following IHowell et al.l (120071) . it may be possible to cali- 
brate out the luminosity difference between the two compo- 
nents with large samples of SNe (e.g., by looking for char- 
acteristic properties in the SN spectra, or by correlating SN 
luminosities with host galaxy types). This will lead to a 
prior constraint on So, although uncertainties will still remain 



on the redshift evolution (as estimated by /d(z)). Assuming 
the two-population fraction, and its redshift evolution, is per- 
fectly known, we analyze the same SN data assuming a Gaus- 
sian prior on 8d- Taking this prior to have zero mean and 
a = (0.25,0.1,0.05), the errors on w are (0.130,0.086,0.072), 
with the best-fit value of w being approximately the same as 
was found for the So = case. Thus, by including a two- 
population systematic in the fit, with a prior on So centered at 
with a 0.05 mag dispersion, we recover the same best-fit w, 
but with ~ 10% degradation in the error bars. 

3.2. Future data 

We now turn to proposed SN surveys, and investigate how 
uncertainties in a possible two-population systematic impact 
measurements of w. We generate mock SN catalogs with 300 
SNe uniformly distributed at z < 0.1, and 2000 SNe in the 
range 0.1 < z < 1.7, similar to a JDEM-like survey dKim et al.l 
120041) . We incorporate an intrinsic Gaussian scatter of 0.1 
mag for each SN, and take the relative fraction of delayed 
and prompt SNe, /d(z), given by Eq. [T|with SB values of A = 
4.4 x 10 _i and B = 2.6. We assume different values (0.025, 
0.05, 0.1 mag) for the underlying two-population bias (So)- 

We fit each mock data set of 2300 SNe, along with the 2 ex- 
isting BAO measurements, to a wCDM model, with the cor- 
responding vvCDM WMAP+HST priors as before. When fit- 
ting the data we consider two extreme cases: 5d = 0, and 5d 
completely unconstrained. Our results for 5d = 0.025, from 
200 separate mocks, are summarized in the right panel of Fig- 
ure 1. The hatched histogram, which peaks at -0.974, depicts 
the case where So = was assumed in the fit. This corre- 
sponds to a systematic being present in the data, but ignored 
in the fit. With an average error on the EOS of 0.029 from 
MCMC, the resulting bias in the best-fit EOS value is ~ \ a 
from the underlying value (w = -1). 

The shaded histogram in the right panel of Figure 1 , having 
a peak at -0.991 with a lcr width of 0.071, shows the dis- 
tribution of the best-fit w with 200 mocks, where we let So 
vary freely while fitting the data to Eq. [3] In this case we 
find no significant bias (w = -1 is recovered within < 0.1c), 
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but the width of the distribution and MCMC results for each 
mock sample show that the errors in the best-fit w increase 
by a factor of ~ 2.5. When the underlying two-population 
bias is So = 0.05 mag, we find w = -0.955 ± 0.030 (assuming 
So = in the fit) and w = -0.986 ± 0.067 (assuming So uncon- 
strained) based on 50 Monte-Carlo realizations. Generating 
data with So = 0.1, and neglecting the systematic in the fit, we 
find w = -0.925 ± 0.029 (50 realizations). As a rough rule, if 
the two-population systematic is neglected, the resulting bias 
in the best-fit w is on the order of the magnitude of the under- 
lying S D . 

Thus far we have assumed no prior knowledge on the val- 
ues of So, although it may be possible to adduce a priori con- 
straints on So through SN population statistics combined with 
correlations to galaxy properties. We model-fit 200 separate 
SN mocks (with intrinsic So = 0.025 mag) with two differ- 
ent priors: S D = 0.025 ± 0.025 and S D = ± 0.025 mag. The 
former case represents knowledge of the true underlying sys- 
tematic (with a la uncertainty of 0.025 mag), while for the 
latter case the central value is incorrectly assumed to be zero 
(with a 0.025 mag uncertainty). With the So prior peaked on 
the correct value (0.025), the distribution of w peaks at -0.993 
with a In width of 0.034 (the Gaussian with dot-dashed line 
in the right panel of Figure 1). With the prior centered on 
the wrong value (0 instead of 0.025), the distribution peaks at 
w = -0.977 (showing a small bias of ~ 0.6er), with the same 
la uncertainty of 0.034. In both cases the errors in w increase 
by ~ 30%, when compared to the equivalent dataset with no 
two-population effect in either the mock data or the fit. 

Thus far we have assumed that we know /d(z) perfectly. 
Even if the two-population model is correct, uncertainties in 
A and B, or equivalently uncertainties in the star-formation 



history, lead to a redshift-dependent uncertainty in fo(z). To 
test the effect of these uncertainties on parameter estimation, 
we generate mock data sets with So = 0.025, and foil) taken 
to be the canonical two-parameter form given in Eq. Q] We 
then fit these data assuming different forms for foil) obtained 
by accounting for uncertainties i n A and B from SB. W e also 
consider an estimate of fo from lAubourg et all (|2007|) . The 
resulting fo curves are shown in the left panel of Figure 2. 
In the right panel we show the resulting bias in dark energy 
EOS measurements, as a result of the uncertainty in fo- The 
uncertainties in the population fraction lead to biases in the 
resulting EOS parameters. To control this bias to the percent 
level, the underlying distribution must be characterized to < 
20%. 

In summary, we have found that a post-calibration shift in 
the standard-candle brightness between delayed and prompt 
SNe can introduce bias in the best-fit dark energy parameters. 
By controlling the magnitude of any resulting two-population 
difference to better than 0.025 mag, the bias can be kept under 
la for a JDEM-liks survey without significantly degrading 
the accuracy of the dark energy measurements. 
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